Luttinger liquid parameter for the spin-1 Heisenberg chain in a magnetic field 
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A recently derived general formula and older numerical results are combined to deduce the be- 
havior of the transverse correlation exponent for the spin-1 Heisenberg antiferromagnetic chain in 
an applied magnetic field: rj ~ 1/2 — (2.0)m -I- 0{m?), where m is the magnetization per site. A 
comparison with the 0(3) non-linear cr-model is also made. 

We consider integer spin Heisenberg antiferromagnetic chains in a magnetic field with Hamiltonian: 

H = Y,[JS^■S,+l~hS'■]. (1) 

i 

In zero field, this model has a singlet groundstate and a gap, A to the lowest excited state which is a spin-triplet of 
magnons. A « .41J for the S = \ Heisenberg model. The gap closes for /i > A at which field the M = 1 magnons 
condense ■1;]. {M is the magnetization.) The magnetization per site, m, is given by the density of these bosons. In 
this condensed phase the correlation functions exhibit power-law decay: 



< s^s;^ > = < sysK > 



(-1) 



3 



The exponent, rj varies continuously with magnetization but was argued, based on earlier results of Haldane 0, to 
behave as: 

r] = 1/2 -0{m). (3) 

The validity of this dilute boson picture of the ground state was confirmed by density matrix renormalization group 
calculations in 0, 0| . In 0] it was shown that the interactions between pairs of dilute magnons could be conveniently 
parameterized by a scattering length, a. This was defined in terms of the phase shift, S{k) for a pair of bosons in 
vacuum, in the small k limit. We write the 2-particle symmetric wave-function as: 

ipix) ^ sm[k\x\ + 6{k)], (|a;| oo), (4) 

where x is the separation of the 2 particles. As A: ^ 0, the phase shift behaves as: 

d{k) -ak, (5) 

defining the scattering length a. (Note that the zero interaction case corresponds to 5{k) = 7r/2 for all fc; k ^ and 
the limit of vanishing interaction strength don't commute.) In some cases a can be determined exactly, as for general 
Bose-Hubbard models In cases like the Heisenberg model where the bosons are collective excitations a, can be 
determined numerically from the finite size corrections to the lowest 2-magnon energy (with either open or periodic 
boundary conditions). With periodic boundary conditions, the lowest 2-magnon excitation energy, (i.e. the energy of 
the lowest state with M = 2) to 0{1/L^) (where L is the system size) is: 

1 /W7r\2 f Aa \ / vn 



E,~-Eo^ 2 VA2 + [7:v/{L - 2aW « 2A + - (— j + j^—j (^—j +0(1/1^). (6) 

[The spin- wave velocity v ~ 2.5 J for the S = 1 Heisenberg model.] A corresponding formula was also derived in 
for the case of open boundary conditions. 

For the S=l chain a was determined by this method to be Q: 

a w -2.0, (7) 

in units of lattice sites, corresponding to about — C/3 where ^ « 6.1 is the correlation length in zero field. [No 
systematic attempt was made to estimate the error in a, which was determined from system sizes L < 100. However, 
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two different estimates, from DMRG data with periodic and open boundary conditions, differed by 7%.] a is the basic 
interaction parameter which determines aU low density properties of the system. In particular, it was shown Q that 
it controls the leading correction, at low density, to the free fermion result for the magnetization density versus field: 



8Aa, 



y/2A{h-A)~^ih-A) 
Sttv 



(8) 



This formula was shown 4] to be in good agreement with numerical data for the above value of a. 

In Q it was shown that a also determines the leading correction to the correlation exponents, at low densities. 
There we expressed the exponent rj in terms of a "Luttinger liquid parameter" g, 

1/^ = 25, (9) 

and derived the result: 

g :^l-2am + 0{m^). (10) 

Combining this exact low density formula for the Luttinger liquid parameter g with the numerically determined 
scattering length leads to the main result of this brief report: 

5 « 1 + 4.0m + O(TO2), (11) 

for the nearest neighbor S — \ Heisenberg antiferromagnet. The corresponding exponents, Eq. might be 

measurable in neutron scattering (or other) experiments. This scattering length is not universal, and is modified 
by crystal field interactions, longer range Heisenberg exchange interactions, etc. On the other hand, our finite size 
method for determining a can be straightforwardly extended to other Hamiltonians. 

This method of determining g at low densities should be much more accurate numerically than conventional tech- 
niques based on direct measurements in the finite ni system. These conventional methods require 

K Af < L/e, (12) 

where L is the system size. On the other hand, a can be determined reliably from studies of the M — 2 sector 
(with periodic boundary conditions, or M = 3 with open boundary conditions) as long as L is large compared to the 
range of the inter-magnon interaction, which is of order ^. (More precisely, we require L ^ 2^ for periodic boundary 
conditions or L ^ 3^ for open boundary conditions.) Numerical results [6j for the Luttinger liquid parameter versus 
magnetization density, g{m), determined from the "Friedel oscillation" exponent using Density Matrix Renormaliza- 
tion Group for a chain of length 120 and 240 are shown in Fig. (1), compared to Eq. 1)11(1 . The agreement at m w .05 
is fair although it seems clear that data for still longer chains would be neccessary to test Eq. (|ll|l . 

The low energy theory for the large-S limit of integer-spin Heisenberg antiferromagnets is given by the 0(3) non- 
linear cr-model (NLctM) (with zero topological term) with Lagrangian density: 

C=^[{l/v){dt$f~v{dJ)% (13) 

with the constraint (f>^ — 1. The exact 2-particle S-matrix for this (integrable) model Q determines 0j the scattering 
length a = — 2^/7r, roughly twice as big as for the S = 1 Heisenberg model. This discrepancy is presumably due to 
the fact that 5 = 1 does not satisfy S ^ 1. Konik and Fendley determined g for all magnetic fields in the NLcrM, 
obtaining in particular a formula for fields slightly above the critical field A: 



25/2 

7r2 VA 



g^l + —^l--l + Oih^A). (14) 



Using the result, a = —2£^/tt, deduced in from the exact S-matrix of 0, together with Eq. ©, and — v/A (which 
is an exact consequence of Lorentz invariance), we see that this result is implied by the general formula for g, in Eq. 

m- 

Finally, we remark that this method of determining the Luttinger liquid parameter at low densities, using Eq. H10() 
and determing the scattering length from studying the 2-boson sector, might also be useful for various other models 
including spin ladders in a magnetic field and bosons in one-dimensional optical lattices, close to the "Tonks gas" 
limit. 
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FIG. 1: Luttinger liquid parameter, g, versus magnetization per site, m. Solid circles are numerical data for a chain of length 
120 and open circles for length 240 from y]; line is Eq. Illll . 
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